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1.  INTRODUCTION 


During  the  design  phase  of  any  new  gun  system,  it  is  important  to  evaluate  the  influence  of  various 
possible  components  and  geometric  configurations.  In  certain  seer  ios,  simulation  of  the  interior  ballistic 
cycle  will  predict  the  growth  of  pressure  wave-  to  amplitudes  w^ach  would  threaten  the  integrity  of  the 
system.  It  is  imperative  to  ensure  the  accuracy  of  these  model  predictions.  Current  interior  ballistic 
simulations  with  Got  .’s  XNOVAKTC  code  (Gough  1990)  (hereafter  called  "XKTC")  can,  for  example, 
predict  propellant  bed  compression  events  where  the  value  of  the  axial  component  of  mixture  stress 
exceeds  100  MPa.  At  the  present  time,  there  are  no  experimental  data  on  mechanical  behavior  of  a 
compacted  bed  of  gim  propellant  grains  at  these  stress  levels,  and  hence  no  way  to  validate  the  model 
predictions.  However,  these  data  should  eventually  be  available  from  a  current  Navy  program  (NSWC/ 
White  Oak  Labs)  which  involves  compression  of  full  scale  (127mm  diameter)  granular  propellant  beds. 

The  present  study  is  an  effort  to  upgrade  the  description  of  mechanical  properties  assigned  to  the 
compacted  propellant  bed  in  the  interior  ballistic  simulation.  This  report  will  attempt  a  brief  review  of 
the  formulation  currently  used  in  the  XKTC  code  (Gough  1990)  and  then  discuss  some  alternative 
formulations  which  might  be  used  to  represent  experimental  data.  Although  an  incompressible  solid  phase 
is  a  basic  assumption  in  the  XKTC  code,  the  discussion  here  will  address  the  implications  of  solid 
propellant  compressibility  on  the  predicted  values  of  wave  propagation  speed  in  a  compacted  aggregate. 

A  number  of  experimental  studies  (Nicolaides,  Wiegand,  and  Pinto  1980;  Lieb  and  Rocchio  1982;  Lieb 
1989;  Gazonas,  Juhasz,  and  Ford  1991;  Gazonas  1991;  Lieb  and  Leadore  1991  &  1992)  have  produced 
valuable  information  concerning  the  mechanical  properties  of  a  single  propellant  grain.  However, 
compression  of  an  aggregate  of  propellant  grains  involves  complex  interactions  at  grain  boundaries  leading 
to  sliding,  deformation,  and  possible  fracture  (depending  upon  the  type  of  material  in  question).  Because 
of  these  interactions,  it  is  not  clear  how  the  mechanical  properties  of  a  single  grain  can  be  used  directly 
to  construct  the  response  of  the  compacted  granular  aggregate.  The  latter  is  required  in  interior  ballistic 
models  which  employ  an  averaging  technique  to  generate  an  elemental  volume  containing  a  binary  mixture 
(i.e.,  that  portion  of  the  volume  not  occupied  by  the  solid  is  assumed  to  be  gas).  An  important 
consequence  of  this  simplification  is  the  necessity  for  a  phenomenological  sub-model  to  link  porosity  of 
the  aggregate  to  its  stress  state.  Oearly,  the  value  of  certain  parameters  used  in  the  description  of  this 
stress  state  must  be  determined  by  matching  the  response  of  the  sub-model  to  data  from  rheological 
experiments  on  compacted  aggregates  composed  of  a  given  solid  material. 
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The  grain- to-grain  "resistance"  forces  associated  with  sliding,  defonnation,  and  fracture  could 
conceivably  introduce  a  strain-rate  dependence  into  the  constitutive  behavior  of  the  aggregate.  At  the 
present  time,  there  is  not  a  great  deal  of  experimental  data  to  guide  a  model  (see  review  by  Conroy  1991). 
Single  grains  of  JA2,  M30  and  XM39  gun  propellant  which  are  subjected  to  compression  in  a  Hopkinson 
Split  Bar  test  (Lieb,  1991)  show  the  presence  of  a  strain-rate  influence,  as  they  do  in  the  Servohydraulic 
tests  (Gazonas  1991)  at  lower  strain  rates.  However,  combustion  of  these  damaged  grains  in  a  closed 
bomb  (Gazonas,  Juhasz,  and  Ford  1991)  indicates  that  the  effective  burning  surface  area  is  nearly 
independent  of  the  strain-rate  history  of  the  sample.  Apparently,  the  only  modeling  attempts  to  account 
for  rate-dependence  within  a  two-phase  mixture  have  been  limited  to  a  "global"  description  (Baer  and 
Nunziato  1986;  Kooker  1990)  which  is  difficult  to  calibrate.  Currently,  all  two-phase  interior  ballistic 
models  known  to  the  authors  have  neglected  rate  dependence  and  assumed  that  changes  in  mixture  density 
are  accommodated  by  an  instantaneous  adjustment  to  the  new  stress  state.  This  will  be  denoted  here  as 
the  "equilibrium"  stress  state.  Note  that  the  values  of  stress  and  associated  wave  propagation  speed 
predicted  by  the  equilibrium  stress  state  are  lower  bounds;  if.  in  fact,  rate  dependence  were  present, 
changes  in  mixture  density  which  occur  at  high-strain  rate  would  lead  to  increased  values.  The  validity 
of  neglecting  rate  dependence  must  be  reevaluated  sometime  in  the  future. 

If  the  propellant  bed  will  be  subjected  to  a  large  amplitude  stress  held,  another  imponant  consideration 
is  the  treatment  of  the  solid-phase  material  itself.  Compaction  of  a  granular  aggregate  can  involve  changes 
in  density  of  the  solid  grains  as  well  as  rearrangement  and  defonnation  of  the  grains.  Quasi-static 
compaction  experiments  (e.g.,  Elban  [1984];  Elban  and  Chiarito  [1986],  and  Sandusky  et  al.  [1988]  arc 
based  on  both  inert  and  energetic  materials  of  small  grain  diameter)  typically  monitor  various  stresses  as 
a  function  of  a  changing  mixture  volume.  To  deduce  the  associated  values  of  porosity  (either  void  volume 
or  solid  volume  fraction)  requires  a  model  of  the  mixture.  If  solid-phase  density  is  assumed  constant 
(incompressible),  then  porosity  values  are  determined  directly  by  the  observed  volume  change.  However, 
for  many  propellants  of  interest  here,  the  incompressibility  assumption  can  lead  to  a  contradiction  when 
analyzing  compaction  data  at  stress  levels  in  the  range  of  100  MPa.  When  the  solid-phase  density  is 
assumed  constant,  the  experimental  values  of  mixture  volume  will  imply  negative  values  of  porosity  (i.e., 
the  mixture  density  has  exceeded  the  theoretical  maximum  density  [TMD]  of  the  solid  at  atmosphenc 
pressure  [denoted  here  as  TMO^]).  Since  accounting  for  solid-phase  compressibility  eliminates  this 
artificial  result,  the  analysis  below  v^ll  examine  the  possible  influence  of  propellant  compressibility. 
Gough’s  (1974)  mixture  theory,  which  has  been  quite  successful  in  simulating  interior  ballistics,  does 
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assume  an  incompressible  solid  phase;  this  is  a  reasonable  assumption  as  long  as  the  magnitude  of  stress 
waves  remains  moderate. 

When  determining  the  coiijiitutive  behavior  of  an  aggregate  of  propellant  grains  with  a  compaction 
experiment,  the  scale  of  experiment  (i.e.,  the  ratio  of  the  test-cell  chamber  diameter  to  the  grain 
diameter)  must  be  consic  ic  J.  A  small-scale  experiment  involving  full-size  propellant  grains  will  create 
a  different  porosity  distribution  and  potentially  a  different  interaction  at  the  circumferential  confining 
bo'mdary.  Since  this  could  introduce  unwanted  distortions  into  the  measured  compaction  response,  it 
would  be  preferable  to  conduct  the  compaction  experiment  at  full  scale.  Although  full-scale  propellant 
bed  compaction  data  is  meager  or  nonexistent  at  the  current  time,  assume  that  such  data  will  be  available 
in  the  future.  If  the  bed  compaction  experiment  is  conducted  under  quasi-static  conditions  (e.g.,  crosshead 
speeds  of  0(10)  cm/min),  then  for  all  practical  purposes,  it  can  be  assumed  that  the  results  have  not  been 
influenced  by  strain  rate.  Thus  the  data  can  be  used  to  construct  or  calibrate  the  rate-independent 
"equilibrium"  stress  state  assigned  to  the  propellant  bed  in  the  interior  ballistic  model  (none  of  the  models 
at  the  present  time  have  a  provision  for  a  strain-rate  effect).  Since  the  NSWC/WO  full-scale  compaction 
experiment  is  to  be  conducted  quasi-statically,  a  reasonable  assumption  is  that  the  propellant  bed  remains 
isothermal.  Furthermore,  the  pressure  of  the  gas  within  the  propellant  bed  can  be  assumed  to  remain  at 
1  atm  since  the  tolerance  between  the  rams  and  the  inside  of  the  cylindrical  chamber  will  not  prevent  the 
escape  of  air  trapped  in  the  chamber. 

2.  MODEL  DESCRIPTION  OF  THE  PROPELLANT  BED 

To  understand  how  the  constitutive  behavior  of  the  compacted  propellant  bed  is  represented  in  the 
XKTC  code  (Gough  1980),  consider  a  two-phase  mixture  composed  of  a  compressible  gas  phase  and  a 
deformable  but  incompressible  solid  phase.  Deformation  of  the  mixture  occurs  in  uniaxial  strain  only, 
along  the  axial  direction,  x.  If  we  assume  here  that  the  mixture  is  confined  within  a  channel  of  uniform 
cross-sectional  area,  then  the  mass  balance  for  the  solid  phase  can  be  written 

{£,).  +  {c.uj,  =  -m/p,  ,  (1) 

where  {  ),  denotes  a  partial  derivative  with  respect  to  x,  and  {  ),  with  respect  to  time.  Although,  in 
general,  the  solid-phase  stress  tensor  will  support  deviator  components,  assume  for  this  discussion  that 


3 


deviator  components  are  negligible:  hence,  the  stress  tensor  can  be  represented  with  the  spherical 
component  (a  pressure),  P,.  The  momentum  balance  for  the  solid  phase  can  be  written 


(2) 


where  D/Dt,  is  the  material  derivative  along  the  solid-phase  streamline.  Now  let  R,  represent  the 
configuration  stress,  or  that  portion  of  P,  greater  than  Pg  due  to  particle-particle  contact  when  the  granular 
solid  phase  is  compacted  into  an  aggregate.  Gough  (1980,  1974)  refers  to  R,  as  the  "intrinsic  average 
intergranular  stress."  With  the  further  constraint  that  the  force  field  at  the  interface  between  phases 
reir.ain.s  in  equilibrium  (i.e.,  P,  =  R,  -i-  Pg),  the  momentum  balance  becomes 


(3) 


Let  o,  be  Gough’s  "nonintrinsic  average  granular  stress"  defined  as 


a,  = 


(4) 


Now  if  R,  can  be  described  by  a  strain-rate-independent  function  of  solid  volume  fraction  (or  porosity), 
then  a,  =  and  it  follows  that 


da, 

de 


i  =  function  (Cj)  3  G(ej)  , 


(5) 


and  hence. 


(6) 


Thus  the  solid-phase  momentum  balance  in  Equation  3  can  be  written 
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(7) 


Du. 


^  g-{e.),  =  f. 


To  visualize  the  role  played  by  the  function  G,  consider  a  compacted,  granular  mixture  of  incompressible 
solids  under  near-vacuum  conditions  (i.e.,  the  contribution  from  the  gas  phase  is  negligible  [Pg  =  0]). 
Then  for  a  small  perturbation  (denoted  by  primes)  about  the  quiescent  state  (denoted  by  subscript  0),  the 
solid-phase  mass  balance  Equation  1  and  momentum  balance  Equation  7  can  be  written,  respectively,  in 
their  linearized  form  as 

(e)  +  c  (u)  =0, 

lO*  ”  ’  ,o\ 


pc  {u}  +  G(c)  =  0. 


Now  differentiating  the  first  with  respect  to  x  and  the  second  with  respect  to  t  leads  to 


which  is  recognized  as  a  wave  equation  with  the  speed  of  propagation  given  by 


(9) 


(10) 


In  the  general,  nonlinear  case,  the  function  G  behaves  as  a  "stiffness”  modulus  which  is  the  product  of 
density  and  propagation  speed  squared. 


G(c.)=P.a.'-  (D) 

In  the  theory  underiying  the  XKTC  code  (Gough  1990),  the  constitutive  assumption  which  defines  the 
solid-phase  stress  tensor  is  embedded  into  the  function  G  (Equation  11)  by  specifying  a  functional 
dependence  of  a,.  The  propagation  speed,  a,,  is  assigned  as 


5 


a,  =  a,  when  loading  , 


(12) 


=  when  unloading  , 


where  the  user-supplied  constant,  aj,  represents  the  speed  of  propagation  during  compressive  loading  when 
the  bed  is  at  the  settling  porosity,  and  the  constant  a2  represents  the  propagation  speed  during 
unloading  from  any  state.  Use  of  Equation  12  in  Equation  1 1  effectively  detennines  the  dependence  of 
a,  on  e,.  To  see  this,  recall  the  definition  of  G  in  Equation  5  which  may  be  written  as 


do. 


=  G(e.)  = 


(13) 


A  straightforward  integration  from  to  (during  loading)  yields  the  result 


So 


e  e 

»  So 


(14) 


which  Gough  refers  to  as  the  "nominal  loading  curve."  Note  that  specifying  a  value  of  a,  uniquely 
determines  o,  as  a  function  of  porosity  of  the  aggregate.  Values  of  stress  predicted  by  the  nominal 
loading  curve  are  sensitive  to  the  value  chosen  for  aj,  particulaiiy  when  the  aggregate  is  compressed  to 
stress  levels  where  porosity  e,  begins  to  vanish.  More  details  of  the  derivation  of  Equation  14  may  be 
found  in  Coiuoy  (1992). 

Wave  structure  in  the  two-phase  mixture  will  be  strongly  influenced  by  the  effective  propagation  speed 
in  the  mixture.  This  wave  propagation  speed  can  be  predicted  from  a  method-of-characteristics  analysis 
of  the  two-phase  equation  system;  the  details  are  beyond  the  scope  of  this  report  but  can  be  found  in 
Gough  (1974).  In  general,  both  solid  and  gas  phases  contribute  to  the  value  of  this  speed.  However,  the 
present  study  is  concerned  with  the  influence  of  the  solid  phase — before  combustion  generates  a  significant 
gas  pressure.  Accordingly,  we  will  suppress  the  contribution  of  the  gas  phase  by  assuming  a  vanishingly 
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small  gas  pressure  (e.g.,  a  compressed  aggregate  under  vacuum  conditions).  Under  these  special 
conditions,  the  effective  speed  of  propagation  in  the  mixture  will  be  a,  given  by  Equation  12. 

3.  COMPRESSIBLE  SOLID  PHASE 


Attempts  to  predict  the  transition  to  detonation  in  granular  energetic  material  have  led  to  the 
development  of  several  two-phase  model  equation  systems  (Baer  and  Nimziato  1986;  Kooker  1990)  which 
account  for  a  compressible  solid  phase.  A  major  difference  with  an  incompressible  model  is  the  addition 
of  a  solid-phase  energy  equation  to  track  a  variable  internal  energy  along  with  the  variable  solid-phase 
density.  Although  the  analysis  by  Kooker  (1990)  has  not  been  used  to  simulate  the  interior  ballistic  cycle, 
this  model  assumes  a  strain-rate-independent  equilibrium  stress  state  for  the  compressed  aggregate  and, 
hence,  should  be  capable  of  predicting  the  influence  of  solid  propellant  compressibility  on  wave  speed. 
Again,  the  scope  of  the  present  study  requires  only  a  brief  summary  of  some  results  from  this  model; 
details  may  be  found  in  Kooker  (1990). 


Similar  to  the  discussion  of  the  Gough  model  above,  assume  here  that  the  solid-phase  stress  tensor 
can  be  represented  with  the  spherical  component  (a  pressure),  P,.  Now,  however,  P,  is  a  function  of 
density  pj  and  specific  internal  energy  e,  (i.e.,  P,(p,,  c,)).  This  function  follows  from  the  assumption  that 
the  Hugoniot  for  the  homogeneous  solid  material  can  be  described  by  a  linear  path  in  the  shock  velocity  - 
particle  velocity  plane,  i.e., 

shock  wave  velocity  =  b^^  +  ♦  particle  velocity  ,  (15) 

which  leads  to 


Ps(Ps-e,)  =  p.  r„ 


r„R 

1  +  “ 

f 

P..V 

r„R 

7  4.  ® 

R  +  1 

R  +  l 

2  J 

R  +  1 

(1-AshRr 


<16) 


where  R  =  (p,  /p.)  -  1  . 


Based  on  this  equation-of-state,  the  bulk  sound  speed  in  the  homogeneous  solid  material  is  then 
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(17) 


When  the  solid  phase  is  assumed  incompressible,  changes  in  porosity  (or  solid  volume  fraction)  of 
the  mixture  are  predicted  directly  by  the  mass  balance  in  Equation  1.  However,  when  the  solid  phase  is 
compressible,  the  mass  balance  becomes 

{  Ps®,  It  +  (  P,E,u,  K  =  -m  ,  (18) 

which  shows  that  changes  in  porosity  and  density  are  coupled.  In  the  Kooker  (1990)  model,  the  additional 
equation  required  to  determine  both  quantities  is  the  constraint  that  the  force  balance  at  the  interface 
between  phases  remains  in  equilibrium  (also  imposed  in  Gough’s  [1974,  1990]  model) 


P.(P.>e.)  -  Pg(Pg.ep  =  p.(c,)  ,  (19) 

where  P,  is  the  configuration  stress  (same  as  R,  in  Gough’s  model)  due  to  particle-particle  contact  when 
the  granular  solid  phase  is  compacted  into  an  aggregate.  Hiis  rate-independent  configuration  stress  is 
assumed  to  be  a  function  of  porosity  as 

P.(e.)  =  -cjl  -  C"'  +  (C”'  -  1)]  In  (l/Cg)  , 

where  C  =  ( 

®  •• 

where  the  values  of  the  parameters  t,,  p,,  P2,  and  B2  are  determined  by  matching  data  from  a  quasi-static 
compaction  experiment  on  the  granular  material  of  interest.  The  momentum  balance  for  the  solid  phase 
will  simply  be  Equation  3  with  R,  replaced  by  P,. 

The  effective  wave  propagation  speed  in  this  mixture  is  also  a  function  of  both  the  gas  phase  and  the 
state  of  stress  of  the  compacted  aggregate.  As  in  the  case  of  Gough’s  (1974, 1990)  model,  the  relationship 
follows  from  a  method-of-characteristics  analysis  of  the  two-phase  equation  system  (Kooker  1990).  Under 
the  same  conditions  postulated  above  for  the  Gough  mode!  (aggregate  in  hydrostatic  stress,  negligible  gas 
pressure),  the  effective  wave  speed  reduces  to 
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a 
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where 


/  'X 

P'’ 

P'-  “a 

V  J 

dP. 
dc,  ’ 


and  a. 
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Ps 

- 

fP.l 
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_ 
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(21) 


In  the  examples  below,  it  is  of  interest  to  compute  the  "compressible"  analog  (denoted  here  as  ajJ  to  the 
incompressible  speed,  a,,  defined  by  (G/pj’^  as  given  in  Equation  11.  Recalling  the  definition  of  G  in 
Equation  5  and  with  o,  =  e,  p,,  then 


G(€,) 

P, 


d(c.P^)/de, 

p; 


(22) 


where  P,  is  defined  by  Equation  20  and  the  values  of  e,  and  must  satisfy  the  equilibrium  constraint 
in  Equation  19. 

4.  DISCUSSION  OF  RESULTS 

The  analysis  in  Section  3  shows  how  the  effective  wave  propagation  speed  in  the  mixture  is  influenced 
by  the  equilibrium  stress  state  of  the  aggregate  of  compacted,  granular  gun  propellant.  Computation  of 
numerical  examples  may  provide  a  better  assessment  of  this  influence.  Since  the  foundation  of  the 
prediction  is  the  equilibrium  stress  state  assigned  to  the  mixture,  the  first  step  should  be  a  comparison 
to  data  from  full-scale,  quasi-static  compaction  experiments  on  various  granular  gun  propellants.  At  the 
present  time,  however,  no  such  data  exist  in  the  stress  range  of  interest  Some  limited  data  at  lower 
amplitudes,  for  example,  have  been  reported  by  Robbins  from  3-in-diameter  cylindrical  compression  cell 
as  described  in  Robbins  and  Conroy  (1992).  Some  of  these  data  for  M30  propellant  at  room  temperature 
are  plotted  in  Figure  1  (Robbins  and  Conroy  1992)  in  the  form  of  applied  force  vs.  a  gas  porosity  which 
is  estimated  from  the  observed  volume  change  and  the  assumption  of  incompressibility.  Figure  1  also 
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includes  theoretical  estimates  of  the  applied  load  which  follow  from  Grough’s  solid-phase  stress  function 
(Equation  14)  for  a  range  of  assumed  values  for  the  aggregate  wave  speed  at  the  settling  porosity. 
Althou^  a  reasonable  estimate  for  the  wave  speed  at  the  settling  porosity  is  a^  =  17S  m/s;  in  this  case. 
Equation  14  does  not  lead  to  a  good  overall  representadcm  of  the  data.  It  is  anticipated  that  full-scale, 
quasi-static  compaction  data  obtained  on  gun  propellants  of  interest  will  likely  be  similar  to  existing  data 
for  small  granulation  material.  Typical  compaction  data  displayed  in  Figure  2  represent  two  very  different 
materials:  Qass  D  HMX  is  a  brittle  ciystalline  explosive  (TMD  of  1.903  g/cm^)  with  an  average  particle 
diameter  of  870  pm,  and  TS-3659  is  a  defonnable  double-base  ball  propellant  (-21.6%  NG  in  NC,  TMD 
of  1.64  g/cm^)  with  neariy  spherical  grains  of  434  pm  diameter.  Details  of  the  experiment  and  the  data 
can  be  found  in  Elban  and  Chiarito  (1986)  for  Gass  D  HMX,  and  in  Sandusky  et  al.  (1988)  for  the  TS- 
3659  ball  propellant 

The  smooth  curves  through  the  quasi-static  compaction  data  shown  in  Figure  2  were  determined  by 
Kooker  (1988,  1990)  as  "best  fits"  based  on  the  assumed  functional  form  for  in  Equation  20  (65.3% 
TMD  HMX:  -  1.0  K  psi,  p^  =  0.,  Bj  -  0.79,  pj  =  1.22;  60%  TMD  TS-3659:  Xj  -  2.6  K  psi,  pj  =  5.5, 
B2  -  3.2,  P2  =  0.08).  Note  that  in  the  early  stage  of  compaction,  the  slope  of  the  mixture  stress  curve 
is  shallow  for  granular  HMX  which  undcigoes  consideraWe  grain  fracture  as  the  aggregate  is  compressed, 
while  the  slope  for  the  deformable  ball  propellant  is  neariy  uniform  at  a  modest  value.  However,  in  the 
later  stage  of  compaction  for  both  materials,  mixture  stress  increases  steeply  as  the  aggregates  are 
compressed  to  densities  near  their  respective  TMD^. 
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Figure  1.  Four  Experimental  Compaction  Curves  (Robbins  and  Conmy  1992'>  for  M3Q  ProneUant  at 
294  K  Compared  to  Gough's  Description  (Equation  U")  for  Values  of  a,  Fgual  to  125.  150. 
175.  and  200  m/s. 


Mixture  Density  (g/cc) 


Figure  2.  Quasi-Static  Compaction  Data  for  Two  Granular  Materials:  65.3%  TMD  Pass  D  HMX. 

870  um  Diameter.  TMD.  =  1.903  g/cm^  (Elban  and  Chiarito  1986).  rb)  59%  TMD  TS-3659 
Double-Base  Ball  Propellant.  434  urn  Diameter.  TMD.  =  1.64  g/cm^  (Sandusky  et  al.  1988'>. 
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If  the  quasi-static  compaction  data  in  Figure  2  represented  actual  gtm  propellant,  then  an  XKTC  (1990) 
simulation  would  first  require  that  the  nominal  loading  curve  of  Equation  14  be  used  to  described  the  data. 
As  an  example  here,  choose  the  HMX  data  which  is  reproduced  in  Figure  3.  Also  shown  in  Figure  3  are 
the  three  curves  which  follow  from  Equation  14  when  aj  is  chosen  to  be  100,  200,  and  300  m/s.  Note 
that  in  the  early  stages  of  compaction,  the  slope  of  the  experimental  data  is  well  approximated  by  the 
curve  based  on  a^  =  100  m/s,  but  in  the  later  stages  the  value  of  an  "effective”  a,  increases.  Although  o, 
given  by  Equation  14  is  not  a  good  approximation  for  the  HMX  data,  it  must  be  remembered  that  Gough 
chose  this  functional  dependence  over  15  years  ago  without  the  benefit  of  any  experimental  data.  From 
that  perspective,  his  choice  was  remarkably  insightful. 


Mixture  Density  (g/cc) 


Figure  3.  Mixture  Stress  vs.  Mixture  Density  for  Isothermal  Quasi-Static  Compaction  of  65.3%  TMD 
Class  D  HMX.  Data  From  Elban  and  Chiarito  (1986):  Solid  Curve  Is  "Best-Fit"  bv  Kooker 
(1988).  Other  Curves  From  Gough’s  Formulation  (Equation  n’t  With  Values  of  From 
100-300  m/s. 
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Predictions  for  mixture  wave  speed  (intentionally  neglecting  the  contribution  from  the  gas  phase)  are 
shown  in  Figure  4,  where  the  solid  line  follows  from  Equation  21  and  the  compressible  analysis.  All  the 
estimates  concur  that  the  effective  propagation  speed  increases  substantially  as  the  porosity  of  the  mixture 
begins  to  vanish  (i.e.,  mixture  density  approaches  that  of  the  TMD  solid).  As  a  result  of  the  fit  to  the 
quasi-static  compaction  data,  Gough’s  incompressible  analysis  (chain-dot  curve  from  Equation  12)  with 
a(  =  100  m/s  increases  slowly  but  then  rises  quite  abmptly  as  mixture  density  aji^roaches  TMD„  of 
1.903  g/cm^.  Retaining  the  same  basic  assumptions  but  representing  the  experimental  data  for  o,  with 
the  solid  line  shown  in  Figure  3  produces  the  "dotted"  curve  in  Figure  4.  Using  the  same  basic  structure 
(i.e..  Equation  22)  but  replacing  p,  and  e,  with  the  "compressible"  values  which  are  roots  of  Equation  19, 
leads  to  the  chain-dash  curve  shown  in  Figure  4.  This  curve  is  a  reasonable  approximation  to  the 
compressible  behavior  up  to  a  mixture  density  of  1.8  g/cm^.  At  mixture  densities  greater  than  1.8  g/cm^ 
in  Figure  4,  the  incompressible  estimates  diverge  significantly  from  the  compressible  prediction. 


Figure  4.  Effective  Wave  Propagation  Speed  in  65.3%  TMD  Qass  D  HMX  During  Isothermal 
Compression  (With  P^  =  0).  Solid  Curve  From  (Eouation  21)  Assuming  Compressible  Solid 
Phase:  Chain-Dot  Curve  From  Goueh  (Eouation  12)  With  a,  =  1(X)  m/s:  Dotted  (Zurve  Is 
Incompressible  fG/o,V^.  Chain-Dash  Curve  Is  Compressible  (G/o,)*^  or  a..  From  Equation  22. 
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Figure  5  gives  a  clearer  picture  of  the  change  in  behavior.  At  low  density,  the  mixture  propagation 
speed  is  controlled  by  the  hindamental  compaction  characteristics  of  the  granular  aggregate  which  are 
closely  approximated  by  a„  (chain-dash  curve  in  Figure  S)  given  by  Equation  22.  At  high  density,  most 
of  the  porosity  has  been  eli  ..:nated  and/or  the  aggregate  has  "locked”  such  that  the  mixture  propagation 
speed  is  dominated  by  the  compressibility  of  the  homogeneous  material  itself  (dotted  curve  in  Figure  S). 
The  compressible  prediction  given  by  Equation  21  transitions  between  these  two  limits. 

The  implication  of  these  results  may  be  easier  to  visualize  by  plotting  mixture  propagation  speed  as 
a  function  of  mixture  stress,  as  shown  in  Figure  6.  For  the  present  example  at  values  of  mixture  stress 
below  SO  MPa,  there  is  little  differeiKe  among  the  predictions  from  the  various  methods.  Important 
differences  appear  near  stress  levels  of  100  MPa.  And  at  2QO  MPa,  there  is  a  dramatic  divergence  such 
that  the  incompressible  result  is  nearly  a  factor  of  two  greater  than  the  compressible  prediction.  In 
general,  as  the  effective  wave  speed  is  reduced,  increases  in  mixture  stress  are  more  likely  to  propagate 
as  shock  waves.  However,  a  more  detailed  investigation  involving  simulations  will  be  required  to  predict 
the  influence  on  the  behavior  in  a  gun  combustion  chamber. 
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Figure  5.  Effective  Wave  Propagation  Soeed  in  65.3%  TMD  Pass  D  HMX  During  Isothermal 
Compression  (With  P^  =  OV  Solid  Curve  From  ^uation  21  Assuming  Compressible  Solid 
Phase:  Chain-Dash  Curve  Is  Compressible  fG/o,l*^  or  a^.  From  Equation  22:  Doned  Curve  Is 
Bulk  Sound  Speed  in  Homogeneous  Solid  (Eouation  171. 
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Figure  6.  Effective  Wave  Propagation  Speed  vs.  Mixture  Stress  in  65.3%  TMD  Pass  D  HMX  During 
Isothermal  Comotession  fWith  =  OV  Solid  Curve  From  Equation  21  Assuming 
Compressible  Solid  Phase:  Chain-Dash  Curve  Is  Compressible  or  a,^  From  Equation 

22:  Chain-Dot  Curve  From  Goueh  fEouation  121  With  a,  =  100  m/s. 

A  similar  analysis  and  comparison  was  made  for  the  double-base  ball  propellant  TS-3659.  The 
experimental  data  from  quasi-static  compaction  shown  in  Figure  2  is  reproduced  in  Figure  7  along  with 
three  curves  which  follow  from  Equation  14  when  a,  is  chosen  to  be  100, 200,  and  300  m/s.  Comparison 
shows  that  the  "effective"  value  of  aj  would  have  to  decrease  considerably  from  its  initial  value  of 
approximately  250  m/s:  the  experimental  data  are  not  well  represented  by  Equation  14.  Again,  this 
illustrates  that  even  a  good  estimate  of  the  compacted  bed  wave  speed  at  the  settling  porosity  does  not 
ensure  that  the  mixture  stress  behavior  has  been  captured. 

Predictions  for  wave  propagation  speed  in  TS-3659  are  displayed  in  Figure  8,  where  again  the  dotted 
curve  follows  fiom  the  incompressible  prediction  similar  to  Gough’s  analysis  but  with  a,  given  by  the 
experimental  data,  the  chain-dashed  curve  follows  from  a^  in  Equation  22  which  incoiporates  the 
compressible  correction  for  p,  and  e,,  and  the  solid  curve  follows  the  compressible  analysis  and 
Equation  21.  The  results  are  similar  to  those  for  the  HMX  case. 
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9%  TMD 


5.  CONCLUSIONS 


All  current  interior  ballistic  calculations  are  based  on  mixture  theories  which  require  a 
phenomenological  submodel  to  relate  porosity  of  the  compacted  aggregate  of  propellant  grains  to  a  stress 
state.  The  accuracy  of  any  simulation  showing  the  growth  of  pressure  waves  to  potentially  dangerous 
amplitudes  will  depend  upon  the  accuracy  of  the  constitutive  behavior  assigned  to  the  compacted 
aggregate.  The  current  formulation  in  the  XKTC  (Gough  1990)  code  assumes  an  incompressible  solid 
phase  and  a  strain-rate  independent  relationship  which  is  uniquely  deflned  by  a  value  of  mixture  wave 
speed  (aj)  at  the  settling  porosity.  Comparison  to  data  from  small-scale,  quasi-static  compression 
experiments  indicates  that  an  improvement  might  be  expected  from  a  direct  representation  of  the  data. 
Estimates  of  the  effective  wave  propagation  speed  in  the  mixture  (with  the  restriction  of  vanishing  gas 
pressure)  show  the  potential  importance  of  including  compressibility  of  the  solid  phase  when  mixture 
stress  approaches  the  level  of  1(X)  MPa.  Accounting  for  solid-phase  compressibility  within  the  current 
interior  ballistic  models  would  require  a  major  restructuring.  However,  a  reasonable  compromise  might 
be  to  alter  the  predicted  incompressible  wave  speed  with  a  correction  for  compressibility. 
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LIST  OF  SYMBOLS 


'sh 


aj 

av 


■'sh 


B, 


•d 

G 


m 


Pi  Slid  p2 


p, 

R. 

TMD 


P, 

Po 

e. 


slope  of  Hugoniot  for  homogeneous  solid  material  (see  Equation  15) 

wave  propagation  speed  in  packed  bed  during  compressive  loading,  when 

wave  propagation  speed  in  packed  bed  during  unloading 

"compressible"  wave  propagation  speed  (see  Equation  21) 

wave  propagation  speed  in  packed  bed  (see  Equation  12) 

wave  propagation  speed  defined  in  Equation  22 

intercept  of  Hugoniot  for  homogeneous  solid  material  (see  Equation  IS) 

parameter  (n-d)  in  expression  for  p,  (see  Equation  20) 

bulk  sound  speed  in  homogeneous  solid  material  (see  Equation  17) 

interphase  drag  force  per  unit  total  volume 

doyde,  as  defined  in  Equation  S 

dCj/de,  evaluated  at  the  settling  porosityn  £,o 

mass  generation  rate  of  gas  (due  to  pyrolysis  or  combustion  of  solid)  per  unit  total  volume 

parameters  (n-d)  in  expression  for  P,  (see  Equation  20) 

static  pressure  in  gas  phase 

spherical  stress  (pressure)  in  solid  phase 

P,  -  Pg  =  configuration  stress  in  Gough’s  analysis 

theoretical  maximum  density  (TMD^  =>  at  atmospheric  pressure) 

velocity  of  the  solid  phase 

configuration  stress  of  packed  bed  =  intrinsic  average  intergranular  stress 
Mie  Gruneisen  coefficient  (assumed  constant  here) 
solid  volume  fraction  =  solid  volumeAotal  volume 
gas  porosity  =  1  -  e,  =  gas  volume/total  volume 
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p, 


X 

{  )x 
{  Ixx 
t 

(  It 
{  la 

e, 

u. 


Dt 


e. 


=  solid-phase  density  =  mass  of  solid/solid  volume 
=  e,  p,  =  non-intrinsic  average  granular  stress 
=  spacial  coordinate 

-  partial  derivative  with  respect  to  space  coordinate 
=  second  partial  derivative  with  respect  to  space  coordinate 
=  time 

=  partial  derivative  with  respect  to  time 
=  second  partial  derivative  with  respect  to  time 
=  small  variance  of  the  solid  volume  fraction  about  a  quiescent  state 
=  small  variance  of  the  velocity  about  a  quiescent  state 
=  settling  solid  volume  fraction 

=  1  -  e,. 

=  1  -  e, 

=  substantial  derivative  along  solid-phase  streamline  =  {  }(  +  u,  { 

=  solid  phase  energy 

=  solid  phase  energy  at  the  settling  porosity 
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Mechanical  Engineering  Department 
ATTN:  0HE20O,  M.  Gerstein 
Los  Angeles,  CA  90089-5199 

1  University  of  Utah 

Department  of  Chemical  Engineering 
ATTN:  A.  Baer 

Salt  Lake  City,  UT  84112-1 194 

1  Washington  State  University 

Department  of  Mechanical  Engineering 
ATTN:  C.T.  Crowe 
Pullman,  WA  99163-5201 
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1  AFELM,  The  Rand  Corporation 

ATTN:  Library  D 
1700  Main  Street 
Santa  Monica.  CA  90401-3297 

1  Arrow  Techtioiogy  Associates,  Inc. 
ATTN:  W.  Hathaway 

P.O.  Box  4218 

South  Burlington,  VT  05401-0042 

3  AAI  Corporation 

ATTN:  J.  Hebert 
J.  Frankie 
D.  Cleveland 
P.O.  Box  126 

Hunt  Valley,  M  D  21 030-0 1 26 

2  Alliant  Techsystems,  Inc. 

ATTN:  R.E.  Tompkins 

J.  Kennedy 
7225  Northland  Dr. 

Brooklyn  Park,  MN  55428 

1  AVCO  Everett  Research  Laboratory 

ATTN:  D.  Stickler 
2385  Revere  Beach  Parkway 
Everett,  MA  02149-5936 

1  General  Applied  Sciences  Lab 

ATTN:  J.  Erdos 
77  Raynor  Ave. 

Ronkonkama,  NY  1 1 779-6649 

1  General  Electric  Company 

Tactical  System  Department 

ATTN:  J.  Mandzy 
1 00  Plastics  Ave. 

Pittsfield,  M A  01201-3698 

1  IITRI 

ATTN:  M.J.  Klein 
10  W.  35th  Street 
Chicago,  IL  60616-3799 

5  Hercules,  Inc. 

Radford  Army  Ammunition  Plant 
ATTN:  D.A.  Worrell 
W.J.  Worrell 
R.  Goft 
C.  Chandler 
L.  Rivenbark 

Radford,  VA  24141-0299 
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2  Hercules,  Inc. 

Allegheny  Ballistics  Laboratory 
ATTN:  William  B.  Walkup 

Thomas  F.  Farabaugh 
P.O.  Box  21 C 

Rocket  Center,  WV  26726 

1  Hercules,  Inc. 

Hercules  Plaza 
ATTN:  B.M.  Riggleman 
Wilmington,  DE  19894 

1  MBR  Research  Inc. 

ATTN:  Dr.  Moshe  Ben-Reuven 
601  Ewing  St.,  Suite  C-22 
Princeton,  NJ  08540 

1  Olin  Corporation 

Badger  Army  Ammunition  Plant 
ATTN:  F.E.  Wolf 
Baraboo,  Wl  53913 

3  Olin  Ordnance 
ATTN:  E.J.  Kirschke 

A.F.  Gonzalez 
D.W.  Worthington 
P.O.  Box  222 

St.  Marks,  FL  32355-0222 

1  Olin  Ordnance 

ATTN:  H.A.  McElroy 
10101  9th  Street,  North 
St.  Petersburg,  FL  33716 

1  Paul  Gough  Associates,  Inc. 

ATTN:  P.S.  Gough 
1048  South  St. 

Portsmouth.  NH  03801-5423 

1  Physics  International  Library 
ATTN:  H.  Wayne  Wampler 
P.O.  Box  5010 

San  Leandro,  CA  94577-0599 

2  Princeton  Combustion  Research 

Laboratories,  Inc. 

ATTN:  N.  Mer 

N.A.  Messina 
Princeton  Corporate  Plaza 
11  Deerpark  Dr.,  Bldg  IV,  Suite  119 
Monmouth  Junction,  NJ  08852 
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3  Rockwell  International 

RocKetdyne  Division 
ATTN;  BA08, 

J.  Flanagan 
J.  Gray 
R.B.  Edelman 
6633  Canoga  Avenue 
Canoga  Park,  CA  91303-2703 

2  Rockwell  International  Science  Center 

ATTN:  Dr.  S.  Chakravarthy 
Dr.  S.  Palaniswamy 
1049  Camino  Dos  Rios 
P.O.  Box  1085 
Thousand  Oaks,  CA  91360 

1  Southwest  Research  Institute 

ATTN;  J.P.  Riegel 
6220  Culebra  Road 
P.O.  Drawer  28510 
San  Antonio,  TX  78228-0510 

1  Sverdrup  Technology,  Inc. 

ATTN:  Dr.  John  Deur 
2001  Aerospace  Parkway 
Brook  Park,  OH  44142 

2  Thiokol  Corporation 
Elkton  Division 
ATTN;  R.  Biddle 

Tech  Library 
P.O.  Box  241 
Elkton,  MD  21921-0241 

1  Veritay  Technology,  Inc. 

ATTN;  E.  Fisher 
4845  Millersport  Hwy. 

East  Amherst,  NY  14501-0305 

1  Universal  Propulsion  Company 

ATTN:  H.J.  McSpadden 
25401  North  Central  Ave. 

Phoenix,  AZ  85027-7837 

1  SRI  International 

Propulsion  Sciences  Division 
ATTN:  Tech  Library 
333  Ravenwood  Avenue 
Menlo  Park,  CA  94025-3493 
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Aberdeen  Proving  Ground 

1  Cdr,  USACSTA 

ATTN:  STECS-PO/R.  Hendricksen 
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1  Emst-Mach-lnstitut 

ATTN;  Dr.  R.  Heiser 
Haupstrasse  18 
Weil  am  Rheim 
Germany 

1  Defence  Research  Agency,  Military 

Division 

ATTN;  C.  Woodley 
RARDE  Fort  Halstead 
Sevenoaks,  Kent,  TNI 4  7BP 
England 

1  School  of  Mechanical,  Materials,  and 

Civil  Engineering 
ATTN;  Dr.  Bryan  Lawton 
Royal  Military  College  of  Science 
Shrivenham,  Swindon,  Wiltshire. 

SN6  SLA 
England 


1  Explosive  Ordnance  Division 

ATTN;  A.  WiWegger-Gaissmaier 
Defence  Science  and  Technology 
Organisation 
P.O.  Box  1750 

Salisbury,  South  Australia  5108 

1  Armaments  Division 

ATTN;  Dr.  J.  Lavigne 
Defence  Research  Establishment 
Valcartier 

2459,  Pie  XI  Blvd.,  North 
P.O.  Box  8800 

Courcelette,  Quebec  GOA  1 RO 
Canada 

1  U.S.  Army  European  Research  Office 

ATTN;  Dr.  Roy  E.  Richenbach 
Box  65 

FPO  New  York  09510-1500 


2  Institut  Saint  Louis 

ATTN;  Dr.  Marc  Giraud 

Dr.  Gunther  Sheets 
Postfach  1260 
7858  Weail  am  Rhein  1 
Germany 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR  80 _ Date  of  Report  February  1993 _ 

2.  Date  Report  Received _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
opernting  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


Ciiy,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


Department  of  the  Army 


OPRCIAL  BUSINESS 


i 


BUSINESS  REPLY  >L\IL 
nRST  cmss  brmit  h  oooi.  apg.  md 

Postage  «ili  &e  oaid  Oy  addressee 


Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-B  (Tech  Lib) 
Aberdeen  Proving  Ground,  MD  21005-5066 


